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Abstract 



We numerically study the statistical properties of differences of zeros of Rie- 
mann zeta function and L-functions predicted by the theory of the efie product. 
In particular, this provides a simple algorithm that computes any non-real Rie- 
mann zeros from very large ones ("self-replicating property of Riemann zeros"). 
Also the algorithm computes the full sequence of non-real zeros of Riemann 
zeta function from the sequence of non-real zeros of any Dirichlet L-function 
("zeros of L-functions know about Riemann zeros"). We also check that the 
first error to the convergence to the classical GUE statistic near is a Fresnel 
distribution. 
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1. Introduction. 

The goal of this article is to check numerically the predictions for the statistics 
of differences of zeros of Riemann zeta function and other L-functions anticipated 
by the theory of the efie ring structure (see [PMlj and [PM2j ). We carry out a 
detailled statistical analysis of the differences of zeros confirming the predicted limit 
distributions. 

The study of the distribution for the differences of zeros is not new. In 1972 H.L. 
Montgomery ( |Molj . |Mo2] ) investigated these statistics in relation with the class 
number problem. Although he states at the beginning of his article jMo2j that 

"Our goal is to investigate the distribution of the differences 'j — 'j' between the 
zeros. " 

he only studies this distribution for differences of nearby (consecutive or semi-locally 
close) zeros. Montgomery formulated the "Pair correlation conjecture" according to 
which, when T — )■ +oo, 

where N{T) ~ ^^TlogT is the asymptotic number of zeros with positive imaginary 
part less than T, < a < P, and the sum runs over the imaginary part of non- 
trivial Riemann zeros. Note the normalizing factor ^ log T that restricts the statistics 
to nearby zeros. According to the well known story, F. Dyson recognized the pair 
correlation distribution for eigenvalues of large random hermitian matrices used by 
physicists (the Gaussian Unitary Ensemble or GUE). This attracted much attention 
since it was believed to give some support to the spectral approach to the Riemann 
Hypothesis proposed by Hilbert and Polya (indeed first by Polya according to |Odj ). 
which asks to identify an hermitian operator whose spectrum is composed by the 
non-real Riemann zeros (this belief is not justified as we will discuss later on). 

Then in the 80 's A.M. Odlyzko performed intensive computations of the distribution 
of normalized differences of consecutive Riemann zeros confirming numerically the 
GUE distribution conjectured by Montgomery ( |0d] ). Later N. Katz and P. Sarnak 
pursued the numerical exploration of zeros for other L-functions ( |KSj ). Their work 
was continued more recently by many others (see jCoj ). 

These authors restrict their study to normalized differences of consecutive zeros. 
Our goal is to study the distribution of global differences at large. These global dif- 
ferences appear on compact sets uniformly distributed in first approximation. At the 
second order they present very significant discrepancies with the uniform distribution 
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(section 12]). The precise location of these discrepancies is highly significant: We no- 
tice a deficiency of differences of non-real Riemann zeros exactly and precisely at the 
same locations as the Riemann zeros. In other words, the statistics of differences of 
Riemann zeros pinpoints the exact location of these zeros. 

In particular, this property indicates that large Riemann zeros do know the location 
of all Riemann zeros, since the statistics remain unchanged by removing a finite 
number (or a density subset) of Riemann zeros. It is indeed checked numerically in 
section [3] that the statistics of very large zeros do find the location the first Riemann 
zeros. 

We also study the differences of the zeros of an arbitrary Dirichlet L-function 
(section H]). We observe the same phenomenon. The differences have a uniform 
distribution except for deficits located exactly at the precise location of Riemann 
zeros. This fact provides a simple algorithm which computes the sequence of Riemann 
zeros from the sequence of zeros of any Dirichlet L-function. This result confirms what 
has been part of the folklore intuition among the community of specialists (see [Coj ): 
Zeros of L-functions do know about Riemann zeros. 

We perform a statistics in sections and E] which is apparently new in the literature. 
We study the differences of the zeros of one Dirichlet L-function Li against the zeros 
of another Dirichlet L-function L2. A similar result is observed: The discrepancies 
from the uniform distribution are located at the zeros of another L-function which 
can be computed explicitly. When L2 is the Riemann zeta function (section [6]), the 
discrepancies are located at the zeros of Li, i.e. Riemann's zeta function plays the 
role of the identity for this mating operation. In section [71 we mate local Euler factors 
with Riemann zeta function. The discrepancies pinpoint to the arithmetic location 
of the poles of the Euler factor. 

We also make a numerical study of the distribution predicted by Montgomery's 
conjecture (section [H]). The theory of the ene product provides a refinement of Mont- 
gomery's conjecture. More precisely, it provides, and we verify numerically, an asymp- 
totics for the error from Montgomery's GUE limit distribution. At the first order we 
observe the predicted Fresnel (or sine cardinal) distribution. These refinements of 
Montgomery's conjecture seem also new. 

The numerical results presented in this article are unexpected without the theory 
of the ene product. Even more surprising is that such elementary statistics have been 
unnoticed so far. We encourage the skeptic readers who want to "put their hand 
upon the wound" to check the numerical results by themselves in their own personal 
computer (any modern personal computer can do the job). It is indeed extremely 
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simple as we indicate below. All statistics presented in this article were performed 
with public software and public data on a low profile regular laptop (of 2005, with 
500 Mb of RAM memory). They are fully reproducible and we provide all necessary 
information to reproduce them. 

Together with the numerical results we give the simple code for the computations. 
These were performed using the public domain statistical software "R". One can 
download and install the program from the site www.r-project.org where tutorials 
are also available. This statistical software runs under Linux and Windows and the 
source code is open. Some source files for zeros of the Riemann zeta function were 
obtained from A. M. Odlyzko's web site |0d] . A large source file for the 35 first 
million non-trivial Riemann zeros and other files with zeros of L-functions are from 
M. Rubinstein public site |Ruj . The author is very grateful to A.M. Odlyzko and M. 
Rubinstein for making available their data. 

We have restricted our statistics to Dirichlet L-functions only because we have only 
access to large lists of zeros for these functions. These results hold in general for more 
general L and zeta functions. We encourage the readers with access to zero data for 
other more general L-functions to perform similar statistics. 

From now on we refer to non-real zeros of Riemann's zeta function simply as "zeros 
of Riemann's zeta function" or "Riemann zeros". The same simplified terminology is 
used for non-real zeros of general L-functions. We also refer to differences as "deltas" . 
We use the notation p for a non-real zero and the notation 7 for p = 1/2+^7. Then the 
sequence of 7's is the sequence of non-trivial zeros of Riemann real analytic function 
on the vertical line = 1/2}. We eventually refer to them also as Riemann zeros. 

Acknowledgments. 

These numerical computations were conducted during 2004 and 2005, and this 
manuscript was written in 2005. It was circulated among some close friends. I am 
very grateful for their support during these years. These results were announced in 
the conference in Honor of the 200th birthday of Galois, " Differential Equations and 
Galois Theory", held at I.H.E.S. in October 2011. 

2. Self- REPLICATING property of Riemann zeros. 

We consider the sequence Z = {pi)i>i of Riemann's zeros with positive imaginary 
part. We label them in increasing order, '^pi > ^pj for i > j. We write pj = 1/2 + i'jj. 
The first few zeros were computed by B. Riemann with high accuracy as was found 
in Riemann's Nachlass ([Ri2]). The tabulation of the first 80 zeros, i.e. all zeros less 
than 200 follows: 
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71 


- 


14.134725142. . . 


72 


= 


21.022039639. . . 


73 


= 


25.010857580... 


74 




30.424876126 . . . 


75 




32.935061588 . . . 


76 




37.586178159 . . . 


77 




40.918719012 . . . 


78 




43.327073281 . . . 


79 




48.005150881 . . . 


7io 




49.773832478 . . . 


7ii 




52.970321478 . . . 


7l2 




56.446247697 . . . 


7l3 




59.347044003 . . . 


7l4 




60.831778525 . . . 


7l5 




65.112544048... 


7l6 




67.079810529 . . . 


7l7 




69.546401711 . . . 


7l8 




72.067157674 . . . 


7l9 




75.704690699 . . . 


720 


= 


77.144840069... 


721 


= 


79.337375020 . . . 


722 


= 


82.910380854... 


723 


— 


84.735492981 . . . 


724 


= 


87.425274613 . . . 


725 


= 


88.809111208... 


726 


= 


92.491899271 . . . 


727 


= 


94.651344041 . . . 


728 


= 


95.870634228 . . . 


729 




98.831194218. . . 


730 




101.317851006.. 


731 




103.725538040 . . 


732 




105.446623052 . . 


733 




107.168611184. . 


734 




111.029535543. . 


735 




111.874659177.. 


736 




114.320220915.. 


737 




116.226680321 . . 
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738 


- 


118 


790782866 . . 


739 


= 


121 


370125002.. 


740 




122 


946829294 . . 


741 




124 


256818554 . . 


742 




127 


516683880 . . 


743 




129 


578704200 . . 


744 




131 


087688531 . . 


745 




133 


497737203 . . 


746 




134 


756509753 . . 


747 




138 


116042055.. 


748 




139 


736208952 . . 


749 




141 


123707404 . . 


750 




143 


111845808.. 


751 




146 


000982487 . . 


752 




147 


422765343 . . 


753 




150 


053520421 . . 


754 




150 


925257612 . . 


755 




153 


024693811 . . 


756 




156 


112909294.. 


757 




157 


597591818.. 


758 


= 


158 


849988171 . . 


759 


= 


161 


188964138.. 


760 


= 


163 


030709687.. 


761 


= 


165 


537069188 . . 


762 


= 


167 


184439978 . . 


763 


= 


169 


094515416. . 


764 


= 


169 


911976479.. 


765 


= 


173 


411536520. . 


766 




174 


754191523.. 


767 




176 


441434298 . . 


768 




178 


377407776 . . 


769 




179 


916484020. . 


770 




182 


207078484 . . 


771 




184 


874467848 . . 


772 




185 


598783678 . . 


773 




187 


228922584 . . 


774 




189 


416158656.. 
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775 = 


192 


026656361 . . . 


776 = 


193 


079726604 . . . 


777 = 


195 


265396680 . . . 


779 = 


196 


876481841 . . . 


780 = 


198 


015309676 . . . 



We use the notation, for 1 < j < i, 

— li ~ Ij ■ 

We fix > 0. Given T > we consider tfie subset of all deltas of zeros: 
AT,iv = {0 < <T-l<j<i<N} c]0, T] 



We denote by Tq = and N — N{Tq). We are interested in the numerical study 
of the distribution of the elements of i^T,N when N — > +oo and T > is kept fixed. 

We represent the histogram of values in At,n by ticks of e = 10~™ for some T > 
and N large. That is, for each integer 1 < A; < 1 + Te"^ we count how many deltas 
5 yield k = [e~^S] (here the brackets denote the integer part). We denote hy Xk > 
this number. The histograms represent the sequence {xk)i<k<i+Te-'^- We represent 
the figures obtained from two statistics. One very fast and the second more intensive 
and precise. The first one, named (a), computes the deltas oi N — 10^ zeros with 
precision e = 10~^ and T — 100. The computation takes about 15 minutes on the 
author's laptop. The statistics named (b) computes the deltas for = 5.10^ with 
precision e = 10^^ and T = 200. This computation took about 2 days on the author's 
laptop. Obviously, reducing T or increasing e greatly diminishes the computation 
time. 



Figure l.a represents the histogram for statistics (a). Figure l.b represents the 
histogram for statistic (b), restricted to the range [0, 100]. Figure l.c represents the 
histogram for statistic (b) in the full range [0, 200]. 
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Deltas of 100 000 Riemann zeros Deltas of 5 million Riemann zeros 




200 400 600 800 1000 2000 4000 6000 8000 10000 

Deltas Deltas 



Deltas of 5 million Riemann zeros 




Figures l.a, l.b and I.e. 



At first sight this distribution seems to converge (once properly normahzed) weakly 
to the uniform distribution on [0,T]. This uniformity for the statistics of global 
deltas of zeros (i.e. non-consecutive) does not seem to have been noted explicitly in 
the literature. It is hinted by the tail asymptotic to density 1 of Montgomery's pair 
correlation distribution, but it doesn't follow from that due to the semi-locality of 
the differences. We discuss this later in more detail. 

But a closer look to the histogram shows some divergence to the uniform distribu- 
tion. The convergence does not appear to be uniform on [0,T]. We can notice one 
major deficiency for small deltas. This appears when we zoom in the picture near 
(see figures 2. a and 2.b for the range of deltas [0, 2]). 
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Deltas of 100 000 Riemann zeros near 



2 4 



10 12 14 16 18 20 
Deltas 



Deltas of 5 million Riemann zeros near 



100 
Deltas 



150 200 



Figures 2. a and 2.b. 



This deficiency is related to the observed fact that the GUE pair correlation distri- 
bution implies that consecutive zeros tend to repel each other. A closer look at the 
figures reveals a compressed scaled GUE pair correlation distribution as expected. 
The factor of compression is 2^1ogTo as one should expect from Montgomery's con- 
jecture. Notice that scaling the figures by the factor 2^1ogTo (as done by those 
authors studying numerically Montgomery's conjectures) pushes away to +00 (when 
To — > -l-oo and N — >■ -|-oo) all the other interesting irregularities of the histogram that 
are the focus of our study. 



Indeed, other divergences to the uniform distribution appear at some special places 
distinct from 0. We notice a remarkable deficit of deltas at certain locations. This can 
be seen clearly by zooming in at several places. Figures 3. a and 3.b shows zooms at the 
interval [10.00,30.00]. Figures 4. a and 4.b arc centered at the interval [30.00,50.00]. 
Figure 5. a and 5.b at the interval [80.00,100.00]. In all these pictures we observe 
at certain precise locations noticeable negative spikes, i.e. a well localized deficit of 
deltas. 
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Deltas of 100 000 Riemann zeros between 10 and 30 Deltas of 5 million Riemann zeros between 10 and 30 




100 120 140 160 180 200 220 240 260 280 300 1000 1500 2000 2500 3000 

Deltas Deltas 



Figures 3. a and 3.b. 



Deltas of 100 000 Riemann zeros between 30 and 50 Deltas of 5 million Riemann zeros between 30 and 50 




300 320 340 360 380 400 420 440 460 480 500 3000 3500 4000 4500 5000 

Deltas Deltas 



Figures 4. a and 4.b. 




Figures 5. a and 5.b. 
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For statistics (b) with T = 200 we can check larger intervals. Figures 6 and 7 are 
centered around the intervals [100.00, 120.00] and [190.00,200.00]. 



Deltas of 5 million RIemann zeros between 100 and 120 Deltas of 5 million RIemann zeros between 190 and 200 




1 0000 1 0500 1 1 000 1 1 500 1 2000 1 9000 1 9200 1 9400 1 9600 1 9800 20000 

Deltas Deltas 



Figures 6 and 7. Statistics (b). 

At this point the reader should take a moment and compare these pictures, and in 
particular the location of the deficiencies, with the tabulated list of Riemann zeros. 

The key observation now is that the location of these negatives spikes is truly 
special. These locations are precisely at the very same location of the Riemann zeros. 
We recognize in figures 3. a and 3.b the locations of the 4 first Riemann zeros. In all 
the Figures 3a, 3b, 4a, 4b, 5a, 5b, 6 and 7 we recognize the location of the zeros in 
the corresponding intervals. Note in particular in figure 6 the two nearby zeros near 
the value 111, 

734 = 111.029535543 ... 735 = 111.874659177 . . . 
We conclude that 

Riemann zeros do repel their deltas. 

This property of the sequence of Riemann zeros is even more surprising considering 
the fact that it is not invariant by translation, i.e. by a global translation of the 
sequence. The set of deltas is independent of such a translation, but obviously not this 
property. The location of each zero is well determined. Any variation on the location 
of a single zero is obviously irrelevant for the distribution of deltas, but the zero will 
then miss the location of the negative spike. Therefore, only the statistics of the deltas 
determines the precise location of the zeros. This implies that any subsequence of 
density 1 of Riemann zeros does determine the whole sequence. For this 
reason we name this property the self repHcating property of the zeros. 
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The self replicating property of the Riemann zeros is completely mysterious without 
the motivation that lies behind this numerical study: The theory of the eiie product. 



We can confirm numerically these observations (in statistics (a) for simplicity) by 
noticing that in the histogram all the deficit values with cumulative count inferior to 
12 500 fall near a Riemann zero, and conversely any Riemann zero yields a group of 
deficit values. The list of the values of k for which < 12 500 in statistics (a) is the 
following: 1, 2, 3, 4, 5, 139, 140, 141, 142, 143, 208, 209, 210, 211, 212, 248, 249, 250, 
251, 252, 302, 303, 304, 305, 306, 327, 328, 329, 330, 374, 375, 376, 377, 407, 408, 
409, 410, 431, 432, 433, 434, 478, 479, 480, 481, 496, 497, 498, 499, 528, 529, 530, 
531, 563, 564, 565, 566, 592, 593, 594, 606, 607, 608, 609, 649, 650, 651, 652, 669, 670, 
671, 694, 695, 696, 718, 719, 720, 721, 722, 755, 756, 757, 758, 770, 771, 772, 792, 793, 
794, 827, 828, 829, 830, 846, 847, 848, 873, 874, 875, 886, 887, 888, 889, 923, 924, 925, 
945, 946, 947, 957, 958, 959, 987, 988, 989. 



To make the main observation more precise, we can, for example, average out 
all deficit values in each group. We discard the first group of values near that 
corresponds to the deficit at (we will come back to this) . Then we find out as many 
groups as Riemann zeros and their averages are denoted by 71,72,73, ■■ ■ They arc 
all very close to the corresponding zero. Table 6 compares the sequence of zeros (7^) 
with the sequence of averages (7,) for all 29 zeros less than 100. We rounded up the 
averages to the first decimal. The matching of the averages 7 with the zeros 7 is 
striking. 
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i 


7'i 


li 


1 


11.131725112 


11.1 


2 


21.022039639 


21.0 


3 


25.010857580 


25.0 


4 


30.424876126 


30.4 


5 


32.935061588 


32.9 


6 


37.586178159 


37.6 


7 


40.918719012 


40.9 


8 


43.327073281 


43.3 


9 


48.005150881 


48.0 


10 


49.773832478 


49.8 


11 


52.970321478 


53.0 


12 


56.446247697 


56.5 


13 


59.347044003 


59.3 


14 


60.831778525 


60.8 


15 


65.112544048 


65.1 


16 


67.079810529 


67.0 


17 


69.546401711 


69.5 


18 


72.067157674 


72.0 


19 


75.704690699 


75.7 


20 


77.144840069 


77.1 


21 


79.337375020 


79.3 


22 


82.910380854 


82.9 


23 


84.735492981 


84.7 


24 


87.425274613 


87.4 


25 


88.809111208 


88.8 


26 


92.491899271 


92.4 


27 


94.651344041 


94.6 


28 


95.870634228 


95.8 


29 


98.831194218 


98.8 



Zeros versus group averages. 



It is also interesting to study the structure of the distribution of the deficit of deltas 
near the zeros. Once properly scaled, we observe a universal distribution for all zeros. 
This distribution is a negative Presnel distribution, i.e. the distribution generated by 
the Presnel integral (also named sine integral, or sine cardinal function) 

smc.(x) = I - 
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We have 



/+00 
sinc7r(x) dx — 1 
■oo 



The Presnel distribution is the Fourier transform 



sinc^(x) = / n(i) , 

of the box function, 

U{x) = H{x + 1/2) - //(x - 1/2) 





for a; 




for \x\ 



We can appreciate this for the histogram plotted in Figure 8. a. Figure 8.b shows a 
more intensive computation with deltas of the first 10 million zeros. 



Deltas of 5 million 



Deltas of 10 million 





Figures 8. a and 8.b. 

In the figures we appreciate a higher frequency noise that blurs the picture. We 
can filter the noise out by standard filtering procedures. The simplest one would be 
to replace (for example) the sequence (xk) by the sequence {fxk) where 



1 



i=-T/2 



The new figures 9. a and 9.b show the pictures with the noise filtered. 
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Filtered deltas of 5 million Filtered deltas of 10 million 




1200 1300 1400 1500 1600 1200 1300 1400 1500 1600 

Deltas Deltas 



Figures 9. a and 9.b. 
Ehe product computation. 

These numerical observations come from the analytic divisor interpretation of the 
eiie product, noted developped in jPMl] and |PM2] . The space of Dirichlet L- 
functions endowed with usual the multiplication and the efie product is a commu- 
tative ring, having a proper normalization of the Riemann zeta function as the ene- 
multiplicative unit. The ene product is associative, not only with respect to multipli- 
cation, but also to infinite arithmetic Euler products. Given two Euler products with 
polynomials (Fp) and (Gp) (with Fp(0) = Gp(0) = 1), 



p 
p 

then 

FiG{s) = l[Fp^Gp{p-') , 

p 

where Fp-kGp is the plain ene product in C of polynomial whose zeros are the product 
of the zeros of Fp with the zeros of Gp, i.e. if 

^.(^) - n (i - f ) 
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then 

a,/3 ^ ^ 

in particular 

(1 - ap-') * (1 - bp-') = 1 - abp-' . 

The main arithmetic property is that for p ^ q,we have logp and log q Q- independent, 
and 

F,ip-') * G,{q-n = 1 . 
Now we denote that for a real analytic function F, 

F{s) = = F{s) . 

The main statistics in this section have its origin in the following computation 

as) * as) = c(s) * as) 

p 

= C(« + l/2)-^ 

R script. 

The following script can be directly feeded into R in order to plot the histogram 
for the deltas for 5 million zeros with precision 10^^ (statistics (b)). The zeros are 
read from the file "zero. data" (one zero per line in increasing order). The reader can 
consult the R tutorial for more elaborate plotting commands. The scripts for the 
other statistics are simple modifications from this one. 

zeroK-scanC'zero.data" , nlines=5000000) 

zero2<-zetal 

x=rep (0,10000) 

N=5000000 

k=0 



for (i in 1:N) 
{ 
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while (((zerol[i]-zero2[i+k])<100.01) & (k+i>l)) 
{ 

k<-k-l 
} 

k=k+l 
j=k 

while ( (zerol[i]-zero2[i+j]>0) & (zerol [i] -zero2 [i+j] <100 . 01) ) 
{ 

d=100* (zerol [i] -zero2 [i+j] ) x [as . integer (d)] =x [as . integer (d)] +1 

} 

} 

barplot (x) 



3. Large Riemann zeros know about all zeros. 

In this section we perform the same statistics as in section 1. but only using deltas 
of large zeros. The convergence is slower, but the results are the same. This indicates 
that zeros with large imaginary part contain full information on the location of all 
zeros. Indeed a density 1 proportion of zeros with large imaginary part contains the 
information about the location of all zeros. 

We perform the statistics with two sets of data. The first one, using Rubinstein's file 
for the first 35 million zeros and selecting one million after the 30th million zero. The 
second and the third are performed with a much larger set of zeros using Odlyzko's 
file containing 10 000 zeros after 10^^ and after 10^^ respectively. The number of zeros 
in these last two statistics is insufficient. These sets of 10 000 zeros are small and 
the distribution of deltas is not even close to the uniform distribution. In the first 
Odlyzko's file all deltas are smaller than 2568 and in the second smaller than 1409. 
Therefore we observe a linear deficit of deltas even for small values of delta when delta 
increases. In order to pinpoint the deficit at the zeros we filter the cumulative data 
on deltas by removing a moving average. Although statistically not as significant as 
the other statistics, the deficit phenomena is still clearly visible at the location of the 
zeros. It is also less visible for large zeros. This indicates a slower convergence. 

The following figures 10a, 10b and 10c are from the first statistics with 10 million 
zeros (7i) with 20.10^ < i < 30.10^. Figures 10b and 10c show the details near the 
top of the uniform distribution. We can appreciate the similarity of these pictures 
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with the previous ones. Figure 10c is centered around the segment [10, 30] and is 
almost identical to figures 3a and 3b. 



Deltas for zeros between 20 and 30 million Deltas for zeros between 20 and 30 million 




100 200 300 400 500 600 700 800 900 50 100 150 200 250 300 350 400 450 500 

Deltas Deltas 



Deltas for zeros between 20 and 30 million 




100 120 140 160 180 200 220 240 260 280 300 
Deltas 



Figures 10. a, lO.b and 10. c. 



In the following figures we illustrate the results for Odlyzko's large zeros near 10^^ 
(statistics (a)) and near 10^^ (statistics (b)). We can observe the linear decreasing 
of the amount of deltas due to the small number of zeros used. We worked with 
Odlyzko's files containing only 10 000 zeros. Paying close attention we can discern 
the deficit of deltas at the location of the zeros. This can be better seen by filtering 
the data by removing a moving average. Figures 12. a and 12. b show that. In Figures 
13. a and 13.b we have the details for deltas smaller than 50. 
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Figures 11. a and 11. b. 
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Figures 12. a and 12. b. 



Filtered deltas for zeros near 10^^12 



Filtered deltas for zeros near 10^^21 




10 15 20 25 30 35 40 45 50 



10 15 20 25 30 35 40 45 50 
Deltas 



Figures 13. a and 13. b. 
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4. Zeros of L-functions do know Riemann zeros. 

In the survey article of B. Conrey on the Riemann Hypothesis ( |Co] ) we can read 
in the section entitled "The conspiracy of L-functions" , 

There is a growing body of evidence that there is a conspiracy among L-functions 
(...) The first clue that zeta- and L-functions even know about each other appears 
perhaps in works of D curing and Heilbronn (...)These results together (...) gave the 
first indication of a connection between the zeros of ({s) and those of L{s,Xd) ■ 

We confirm in this section that zeros of L-functions do know about all Riemann 
zeros. Indeed we provide a simple algorithm that builds the sequence of Riemann 
zeros from the sequence of zeros of any Dirichlet L-functions. Our first example is 
for the simplest non-trivial L-function: We show how to recover Riemann zeros from 
the zeros of L^.^, where xs is the only character of conductor 3. 

We perform the statistics for the deltas of the zeros of L^.^ as done in section [2] 
for Riemann zeros. This time we observe that the deficit values for the deltas of 
zeros of L^g is located precisely at Riemann zeros. As in section [2] we perform one 
statistic with 100 000 zeros of L^.^ and precision 0.1 for the deltas, and another, more 
intensive, with 5 million zeros of L^.^. Figures 14 show the histogram of deltas for 
both statistics. Figures 15 show the details in the interval [10,30], figures 16 for 
[30, 50], and figures 17 for [80, 100]. We observe in figures 18 the deficit of deltas near 
verifying Montgomery's prediction. 

The similarity of these figures with those in section [2] is clear. Recall though that 
they are generated from a very different set of data. 




Figures 14. a and 14. b. 
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Deltas of 100 000 L_3 zeros between 10 and 30 



Deltas of 1 million zeros of L 3 between 10 and 30 
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Figures 15. a and 15. b. 



Deltas of 100 000 L 3 zeros between 30 and 50 



Deltas of 1 million zeros of L_3 between 30 and 50 
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Figures 16. a and 16. b. 



Deltas of 100 000 L_3 zeros between 80 and 100 



Deltas of 1 million zeros of L 3 between 80 and 100 
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Figures 17.a and 17.b. 
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Deltas of 100 000 L_3 zeros near 



Deltas of 1 million zeros of L 3 near 
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Figures 18. a and 18. b. 



Now we perform the same statistics for the zeros of other Dirichlet L-functions L^. 
We perform the statistics for the deltas of 1 miUion zeros, for deltas in [0, 100], and 
with precision 0.1. This time we consider a real and a complex non-real character. For 
a complex non-real character, the associated L-function is not real analytic, and the 
zeros are no longer symmetric with respect to the real axes. Therefore we compute 
the deltas of those with positive imaginary part and the deltas for those with negative 
imaginary part, and we compute the cumulative result. Since the sequence of zeros 
is not symmetric with respect to 0, we take the first million zeros in the following 
sense: We order the zeros by absolute value and we consider the first million of them 
for the statistics. 

The first statistics is for x = X4: the only primitive character of conductor 4. The 
character is real and the associated Dirichlet function real- analytic. The second 
statistics is for x = X7,3) o^ie of the primitive complex characters of conductor 7. 
Figures 19 show the histograms of deltas in [10, 30]. Figures 20 show the histograms 
of deltas in [30,50]. Again we find that the deficit locations coincide with Riemann 
zeros. 
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Deltas of 1 million zeros of L 4 




Deltas of 1 million zeros of L 7 
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Figures 19. a and 19. b. 




Figures 20. a and 20. b. 



We analyze next the distribution of deltas near 0. We plot the histogram near of 
the deltas of 1 million zeros with precision 0.01. We observe the predicted GUE pair 
correlation distribution as pictures 21 show. 



24 



R. PEREZ MARCO 



Deltas of 1 million zeros of L 4 near 



Deltas of 1 million zeros of L 7 near 
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Figures 21. a and 21. b. 



Eihe product explanation. 

We note that for a real character x, 



The ene product explanation of the first numerical result for X3 is based on the 
following computation (where we denote by Xo the principal character modulo 3) 



^X3 * ^X3 ~ ^X3 * ^X3 



=n(i-^3(p)wi 

= 11(1 - Xo{p)p-'^'-') 
p 

^(l_3-i/2-3)-i-Q(i_^-i/2-.) 
p 

^(l_3-i/2-.)-i^(g + l/2)-^ . 

In general, for an arbitrary character x modulo n, we recognize the distribution of 
the deltas of the zeros of in the result of the eiie product of with L^. We have 
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Also observe that 

X-X = Ixl^ = Xo , 

where xo is the principal character modulo n. 
Therefore we have 

-\\{{.^-x{p)p")A * 1 11(1 -^('^)'^"'^)^' 

= - xWT^ * (1 - xwr^ 

V 

=n(i-ix(p)iv^vi 

V 
V 

Observe that the zeros of each Euler factor 
are for k ^1^, 

1 . 27r , 

2 logp 

According to the explanation with the eiie product, we should observe a deficit of 
deltas (with a lower order amplitude) near locations multiples of the fundamental 
harmonic 



logp 

for A; e Z. 

For p — 2, p — 2) and p — 7 we have, 
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= 9.0647... 



log 2 

-^ = 5.7192... 

log 3 

3.2289 . . . , 



log 7 

A good eye can spot a trace of these deficit locations in the figures 16, 17, 19 and 
20. In particular comparing these figures with figures 3 and 4. 

Script. 

We provide the script for computing the deltas of the zeros of a complex non-real 
L-function since it is slightly different from the previous ones. Here we feed the 
program by reading into Rubinstein's table "zeros-0007-2000000" which contains the 
first zeros for each of the three primitive characters of conductor 7. The zeros for the 
complex character that we are considering are those after row 2 000 000. 

zeros<-read . table ( "zeros-0007-2000000" , skip=2000000 ,nrows=1000000) 

z<-zeros [,3] 

z.plus<-z [z>0] 

z . minus< — z [z<0] 

x=rep (0,10000) 

zetal<-z .plus 

zeta2<-z .plus 

N=length(z .plus) 

k=0 

for (i in 1:N) 
{ 

while ( ((zetal[i]-zeta2[i+k])<100.01) & (k+i >1) ) 
{ 

k<-k-l 
} 

k=k+l 
j=k 

while ( (zetal [i]-zeta2[i+j]>0) & (zetal [i] -zeta2 [i+j] <5 . 01) ) 
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{ 

d=100* (zetal [i] -zeta2 [i+j] ) 

X [as . integer (d)] =x [as . integer (d)] +1 

} 

} 

zetal=numeric() 
zeta2=numeric() 
zetal<-z . minus 
zeta2<-z .minus 
N=length(z .minus) 
k=0 

for (i in 1:N) 
{ 

while ( ( (zetal [i]-zeta2[i+k])<100. 01) & (k+i >1) ) 
{ 

k<-k-l 

} 

k=k+l 
j=k 

while ( (zetal [i]-zeta2 [i+j] >0) & (zetal [i] -zeta2 [i+j] <5 . 01) ) 

{ 

d=100*(zeta[i] -zeta3[i+j] ) 

X [as . integer (d)] =x [as . integer (d)] +1 

} 
} 

5. Zeros of L-functions replicate from their mating with Riemann 

ZEROS. 

We present in this section and the next one a new type of statistics. We do study 
the statistics of differences of zeros of an L-function L^j with the zeros of another 
L- function L^^. We name this operation the "mating" of zeros of L^^ and L^^. As 
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predicted by the ene product theory, it appears that the sequence of Riemann zeros 
plays the role of the unit for this mating operation. More precisely, the statistics of 
this section verify that the mating of Riemann zeros with those of another L-function 
L yield as deficit values the zeros of L itself. 

We perform the statistics mating the Riemann zeros with the zeros of L^,^ where 
X3 is as before the only primitive character of conductor 3. The function L^g is 
real analytic and its zeros are symmetric with respect to the real axes. We consider 
only the non-real (i.e. non-trivial) zeros with positive imaginary part. We denote by 
(1/2 + i7-^^)i>i, or simply (7i"'''*)j>i, the zeros of Ly,^, with i i->- 7^-^^ increasing.. The 
first 18 ones, less than 51, are the following 



(3) 

7i 




8.039737156 


(3) 
72 




11 


.24920621 


(3) 
73 




15, 


.70461918 


(3) 

7a 




18. 


,2619975 . 


(3) 
75 




20. 


,45577081 


(3) 
76 




24. 


,05941486 


(3) 
77 




26. 


,57786874 


(3) 
78 




28. 


,21816451 


(3) 
79 




30. 


,74504026 


(3) 

7io 




33. 


,89738893 


(3) 

7ii 




35. 


,60841265 


(3) 
7i2 




37. 


,55179656 


(3) 
7i3 




39. 


,48520726 


(3) 

Tu 




42. 


,61637923 


(3) 
7i5 




44. 


,12057291 


(3) 
7i6 




46. 


,27411802 


(3) 
7i7 




47. 


,51410451 


(3) 
7i8 




50. 


,37513865 
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This time the "deltas" are differences 



^Lj = 7i - 7 



(3) 



We perform statistics (a) with 1 <i,j < 100 000 and statistics (b) with 1 <i,j < 
10^. We look at deltas in [0,50] with precision 0.1. The results are presented in the 
following figures. 



Mating of 100000 Riemann and L_3 zeros 



Mating of one million Riemann and L_3 zeros 
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Figures 22. a and 22. b. 



Mating of 100000 Riemann and L_3 zeros 



■Mating of one miiiion Riemann and L_3 zeros 
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Figures 23. a and 23. b. 
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Mating of 100000 Riemann and L_3 zeros Mating of one miiiion Riemann and L_3 zeros 
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Figures 24.a and 24.b. 



We observe that this time the deficient locations for the deltas happen exactly at 
the location of the zeros of L^.^. We easily recognize in figures 23. a and 23. b the 
location of the first zeros of L^g. We can check the full list of zeros less than 50 by 
looking also at the figures 24. a and 24. b. We conclude that the zeros of L-functions 
replicate mating them with Riemann zeros. 

A new feature is that near we no longer have a GUE distribution for the deltas. 
As the theory of the eiie product explains, the deficit at only occurs when we have 
symmetric zeros, i.e. we mate the zeros of L^j with those of L^^ when 



Xi = X2 , 



and we have an atomic mass at that comes from the sum of symmetric zeros of L^j 
and L^^. Thus if the character is not real, then we don't have a GUE distribution, 
not even a deficit, but the Riemann Hypothesis is still conjectured, thus there is no 
direct relation between the Riemann Hypothesis and Montgomery Gonjecture. The 
author knows no reference in the literature for this observation. 
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Mating of 100000 Riemann and L_3 zeros 



Mating of one miiiion Riemann and L_3 zeros 



Figures 25. a and 25. b. 



Next we perform the same mating statistics of Riemann zeros with zeros of L^^^. 
Recall that the zeros of this non-real analytic L-functions arc not symmetric with 
respect to 0. We perform two statistics. We consider the first 100 000 Riemann zeros 
and compute all deltas with positive (resp. negative taking their negative value) zeros 
of L^^ 3 . The list of the first positive zeros of L^^ 3 less than 50 are 
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(7+) 

7i 


= 4.356402 . 
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The list of the first negative zeros 



of Lp^^ 3 less than 51 are 
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(7-) 

7i 
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The following figures show the result of the numerical statistics. Figures (a), resp. 
(b), are for the mating against positive, resp. negative, zeros. We observe for sta- 
tistics (a) that the deficient locations do correspond to values of the positive zeros. 
For statistics (b) we observe these locations at the values of the negative zeros. In 
particular in figures 27 we appreciate the location of the first zeros. 




Figures 26. a and 26. b. 




Figures 27.a and 27.b. 
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Mating lO'^S Riemann zeros with positive L_7 zeros Mating 10'^5 Riemann zeros witli negative L_7 zeros 




Figures 28. a and 28. b. 

We perform a final statistic in order to check the nonexistence of the GUE distri- 
bution, and not even a deficit of deltas effect near 0. We compute the deltas near 
against half million Riemann zeros with precision 0.01. The results are shown in 
Figures 29, figure 29. a for positive deltas and figure 29. b for negative ones. Figures 
29 show the deltas with double precision in the range [0, 2]. The reader can compare 
directly these figures with Figures 21. c and 21. d. The conclusion is clear: No GUE 
distribution near 0. 



Mating of 5.10'^5 Riemann zeros with positive L_7 zeros iVIating of S.IO'^S Riemann zeros witli negative L_7 zeros 
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Figures 29. a and 29. b. 

We check also from these statistics the location with double precision the first 
positive and negative zero in figures 30. a and 30. b which show the deltas in the range 
[0, 10] with double precision. We can appreciate distinctly with double precision in 
figure 30. a both positive zeros less than 10 

7^ +) = 4.356402 . . . 7^ +^ = 8.785555 . . . 
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and in figure 30. b both "negative" zeros less than 10 

7^ = 6.201230 . . . 7^ = 7.927431 . . . 



Mating of 5.10'^5 Riemann zeros with positive L_7 zeros KAating of 5.10^5 Riemann zeros with negative L_7 zeros 
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Figures 30. a and 30. b. 
Ene product explanation. 

The computation follows. We have for any Dirichlet L- function L^, 

-U(^-x{p)p-'^'p-n 
p 

^ll(^-x{p)p-^''''^'^) 

p 

^L^{s + 1/2)-' . 

Therefore we recognize that the mating of zeros of with Riemann zeros has 
deficient deltas at the location corresponding to the imaginary part of zeros of L^. 

Scripts. 

Below is the script we used in order to produce the previous figures. The script is 
slightly different from previous ones since we compute separately positive and negative 
zeros. The cumulative positive deltas are stored in the list "x" and the negative in the 
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list "y". The zeros of Dirichlet L- function of conductor 7 are stored in Rubinstein's 
file "zeros-0007-2000000", and Riemann zeros are from Rubinstein's file "zeros-0001- 
35161820". 

zerosL<-read . table ( "zeros-0007-2000000" , skip=2000000 , nrows=1000000) 
z<-zerosL [,3] 
z .plus<-z [z>0] 
z . ininus< — z [z<0] 

zerosR<-scan("zeros-0001-35161820" ,skip=0,nlines=100000) 

x=rep (0,500) 

zetaK-zerosR 

zeta2<-z . minus 

N=length(zeta) j=l 

for (i in 1:N) 

{ 

while ( (zetal[i]-zeta2[j])>50.1 ) 
{ 

} 

1=0 

while ( ((zetal[i]-zeta2[j+l])>0) & ((zetal [i]-zeta2[j+l] )<50. 1) ) 
{ 

d=10* (zeta [i] -zeta3 [j +1] ) 

X [as . integer (d) ] =x [as . integer (d) ] +1 

1=1+1 

} 
} 

y=rep (0,500) 
zetaK-zerosR 
zeta3<-z .plus 
N=length(zeta) 

for (i in 1:N) 
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{ 

while ( (zetal [i]-zeta3[j] )>50. 1 ) j=j+l 
1=0 

while ( ( (zetal [i]-zeta3[j+l])>0) & ((zetal [i]-zeta3[j+l] )<50. 1) ) 
{ 

d=10* (zetal [1] -zetaS [j+1] ) 

y [as . integer (d)] =y [as . integer (d)] +1 

1=1+1 

} 
} 



6. Mating of general L-functions. 

In this section we perform similar statistics to those in the previous section but 
mating the zeros of two Dirichlet L-functions and L^^. This time we observe 
that the deficient locations for the statistics of deltas correspond to the zeros of an 
arithmetically well determined, namely L^^^^. 

For a character x we denote by its conductor. We have 

fx~ fx ■ 

All characters considered are primitive, i.e. defined modulo its conductor. Let Xi ^^'^ 
Xi be two characters. If /^^ A — 1 then the conductor of X\Xi is 

/xiX2 ~ /xi-/x2 • 

The first complex non-real Dirichlet character has conductor 5. Therefore the mat- 
ing of two Dirichlet L-functions of complex non-real characters with distinct conduc- 
tors has conductor at least 35. We have only access to Rubinstein's public data that 
contains large files of zeros for Dirichlet L-functions with conductor < 19. Therefore 
we limit our numerical computation to real characters for which we can check the 
result with the available data. This is done only for checking purposes. Note that we 
could indeed compute, with a rough precision, the zeros of higher conductor Dirichlet 
L-functions (for example 35) by using Rubinstein's data of conductors < 19. 

We choose to mate the zeros of Dirichlet L-functions L^g of conductor 3, and L^^ of 
conductor 4. We should obtain the zeros of the only Dirichlet L-function of conductor 
12, L^j2- The hst of the first zeros of L^^^ less than 50 is 
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Mating S.IO'^S zeros of L_3 with L_4 
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Figure 31. 



Mating S.IO'^S zeros of L_3 with L_4 
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Figure 32. 



Mating 5.1 0*^5 zeros of L 3 with L 4 
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Figure 33. 
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Mating S-IC^S zeros of L_3 with L_4 
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Figure 34. 



Again in this situation there is no GUE distribution near since the zeros of L^g 
and L^^ are not symmetric. Figure 35 shows the histogram of the deltas in the range 
[0, 2] with precision 0.01. This figure is to be compared to figures 21. 



Mating 5.10^5 zeros of L_3 with L_4 
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Figure 35. 



Ene product explanation. 



The computation follows. We have for any pair of Dirichlet L-function L^^ and 
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TT(i-xi(p)p"')"' 

P 


* TT(i - X2(g)r' 


p 


*(1-X2(P)P-^)-' 


Yli^ - {XiX2){p)p~ 
p 


1/2^-) 


Yli^ - iXiX2){p)p~ 




p 





s\-l 



7. Mating with local Euler factors. 

In this section we study the mating of L-functions with local Euler factors. We 
mate the positive imaginary part of Riemann zeros with the sequence of positive 
imaginary part of zeros of Euler factor 

/,(.) = (l-p-) , 

which is the arithmetic sequence 

Ap) _ 27r 
"logb)""' 

with A; G Z. 

More generally we can consider the mating with general Euler Dirichlet local factors 

f,Js) = (1 - xip)p-') , 

but we restrict the numerical statistics to Riemann Euler local factors. 

In order to have a substantial number of deltas we need to work with a very large 
file of Riemann zeros because this time the arithmetical sequence of Euler zeros has 
constant density. For a given number of Riemann zeros we pick more deltas if the 
prime p is large. For this reason, with our limited data of zeros available, we run the 
statistics for p — 23, 

271 

— -— = 2.00389 . . . , 

log(23) 

and the first 10 million Riemann zeros (statistics (a)); and for p = 67, 

27r 

-—— = 1.494327... , 
log(67) 

and the first 15 miUion of Riemann zeros (statistics (b)). 
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The histograms obtained (figures below for different ranges of deltas) show a uni- 
form distribution with deficit locations at the corresponding Euler zeros, for /c G Z, 

The deficit of deltas at these locations can also be seen for statistics with smaller 
values of p, but because of the "small" number of Riemann zeros at our disposition, 
the statistics is too poor 



Mating of Riemann zeta with Euler factor p=23 



Mating of zeta with Euler factor p=67 



.11 III, III III. II 
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Figure 36. a and 36. b. 




Figure 3 7. a and 37.b. 



Also our laptop has a RAM memory of 512 Mb which does not allow R to read a vector with 
more than 20 million Riemann zeros. 
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Figure 38. a and 38. b. 



Ene product explanation. 

We have 



g 

= (l-p-^)*(l-p-^) 

= fp{s + 1/2) . 
We have in general for an arbitrary Dirichlet L-function: 

^ {I - x{p)p-T' ^ - p'T' 
= (1 - xip)p-'^'p-') 

= f,Js + 1/2) . 
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8. Fine structure of deltas near 0. 

We come back in this section to analyze the fine structure of the statistics of deltas 
near 0. 

As we have observed, Montgomery conjecture is verified numerically for the deltas 
of zeros of arbitrary L- functions, but not for the mating of zeros of non-conjugate 
L-functions. As explained, the GUE distribution arises because of the symmetry of 
the zeros mated, i.e. it is a genuine real-analytic phenomenon. 

For the Ricmann zeta function, the ciic product analysis of the distribution of 
deltas near zero reveals that after the first order GUE correction, we have a second 
order term corresponding to the pole of C{s + 1/2) at s = 0. This yields a positive 
Presnel distribution. 

We verify numerically this first order correction to the GUE distribution. In the 
numerical application we consider the deltas of 5 million Riemann zeros. These have 
imaginary part less than Tq for To — 2 630 122. We consider the frequency 

cuo = ^ log To = 0.02352714 . . . 
We correct the histogram of the deltas by adding a GUE density 



t ^ A 



/ sin(7ra;oi) Y 



Figure 39 shows the distribution of the deltas and figure 40 the corrected distribu- 
tion of the deltas showing the Presnel distribution (both in the range [0, 2]). Figures 

41 and 42 show the same data in the range [0, 6] for a better view of the queue. We 
have adjusted the coefficient A in order to fit the best with a Fresnel distribution 
looking at the first minima and the second maximum. 



Deltas of 5 million Riemann zeros near 
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Figure 39. 
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First correction deita distribution near 
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Figure 40. 



First correction delta distribution near 
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Figure 41. 



Deltas of 5 million Riemann zeros near 
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Figure 42. 
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Script. 

We have previously stored tiie result of computing the deltas of 5 million Riemann 
zeros in the range [0, 200] with precision 0.01. The variable x contains the cumulative 
count of deltas with the stated precision. Note also that tt.wo/IOO = 0.0739127 . . . 

loadC'data 5 million") 

correction=numeric() 

for (k in 1:200000) correction[k]=139000*(sin(0.0739127*k)/(0.0739127*k))2 

xK-x+correction 

barplot(xl) 
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